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Abstract 

Dynamical spontaneous breaking of some discrete symmetries including spe- 
cial parities and time reversal and their restoration at finite temperature T are 
researched in 3D Gross-Neveu model by means of Schwinger-Dyson equation in 
the real-time thermal field theory in the fermion bubble diagram approximation. 
When the momentum cut-off A is large enough, the equation of critical chemical 
potential fj, c and critical temperature T c will be A-independent and identical to 
the one obtained by auxialiary scalar field approach. The dynamical fermion mass 
m = m(T, //), as the order parameter of symmetry breaking, has the same (T c — T) 1 / 2 
behavior at T ~ T c as one in 4D NJL-model and this shows the second-order phase 
transition feature of the symmetry restoration at T > T c . It is also proven that no 
scalar bound state could exist in this model. 
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The Gross-Neveu (GN) model [1] in 2+1 space-time dimension is a model with four- 
fermion interactions. As a well-known and applicable model to codensed matter physics, 
it has been extensively researched [2-6] . The coventional approach is to introduce an aux- 
ialiary scalar field, then to replace the four-fermion interactions by the Yukawa couplings 
between the auxialiary scalar field and the fermion fields. Alternatively, one can also use 
Schwinger-Dyson equation directly to research the fermion condensates induced by the 
four-fermion interactions and relevant dynamical behavier. In the latter approach, no ef- 
fective scalar field is introduced and, as far as some problems, for example, the calculation 
of the propagators for the bound states composed of fermions at finite temperature, are 
concerned, this approach could be more straightforward and convenient. On the other 
hand, almost all the discussions about the finite temperature behavior of the model were 
conducted in the imaginary-time formalism of thermal field theory. However, we also 
want to examine the application of the real-time formalism of thermal field theory to 
the model. Therefore, in this letter, instead of using the auxialiary scalar field and the 
imaginary-time thermal field theory, we will research the 3D GN model at finite temper- 
ature by means of Schwinger-Dyson equation in the real-time thermal field theory. All 
the discussions will be made in the fermion bubble diagram approximation. This approx- 
imation amounts to the leading order oil/N expansion and is sometimes also called the 
mean field approximation. Since 3D GN model is renormalizable in 1/N expansion [5] 
and in 2+1 dimension, there is no the kink effect which could appear in a 1+1 dimension 
system [7], the mean field approximation can be considered as a reliable one. The results 
will also be compared with the previous works' ones based on the auxialiary field and the 
imaginary-time thermal field theory. 

The Lagrangian of the model can be written by 

TV TV 

c(x) = j;^#) + ^E^( a; )^W] 2 - (!) 
k=i ziv k=i 

In 3 dimensions, the coupling constant g has mass dimension -1 thus the model is pertur- 
batively non-renormalizable. However, this fact does not affect our discussions based on 
the 1/N expansion. The 7 M (fx = 0, 1,2) can be taken as 2 x 2 matrices and related to 
the Pauli matries a l {i = 1, 2, 3) by 

which satisfy the Dirac algebra 

{ 7 M 7 "} = 2«r, ^ = 0,1,2. (3) 

In this representation ipk{x) is two-component complex spinor with N "color" degrees of 
freedom. Notice that there is no 2 x 2 matrix which anti-commutes with all the three 7 M , 
i.e. no "75" matrix exists in 3 dimensions. Consequently C(x) has no chiral symmetry, 
though ipk{x) are zero- mass fields. However, it is easy to verify that the action of system 
is invariant for the discrete transformation (omitting the index k of ip) 

^(x)^^(x), if^) = (^g), (4) 
the special V\ and Vi parity transformations 

^t,x\x 2 )^^(t,-x\x 2 ), (5) 

^t,x\x 2 )^^(t,x\-x 2 ), (6) 
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and the time reversal 

VKM\* 2 )^7>(-M\* 2 )- (7) 

We indicate that the mass term —mip(x)ip(x) (if it exists) is not invariant under Z 2 , Vi, 
V 2 and T. Hence, if the four-fermion interactions could lead to generation of a non-zero 
dynamical mass of the ip- fields, then the above Z 2 , Pi, V 2 and T discrete symmetries will 
be spontaneously broken. Since the broken symmetries are discrete, no Goldstone bosons 
are expected to appear. In the following we will research generation of the dynamical 
mass of ip{x) fields and its finite temperature behavior in the fermionic bubble graph ap- 
proximation. This approximation, as the leading order of the 1/N expansion, will be able 
to give us a clear insight into the breaking and restoration of above descrete symmetries. 

At zero-temperature T = 0, assume the four-fermion scalar interactions g J2k=i (ffiipk) /2A/" 
can lead to the fermion condensates X^Li^^fe) 7^ 0, then we will get the dynamical 
fermion mass m(0) = —(g/N) X^Li^^fc) which satisfies the gap equation 



idH 

{2nf I 2 - m 2 (0) + ie 



19 I , n , 3 n ■ ~ (8) 



After Wick rotation, angular integration and introduction of 3D Euclidean momentum 
cut-off A, Eq. (8) becomes 



g_A 

TT 2 



, m(0) A 
1 - — arctan ■ 



A m(0) 



(9) 



We indicate that since A corredponds to the square root of the 3d squared Euclidean 
momentum, its introduction does not explicitly break the discussed discrete symmetries 
including Vi, V 2 and T. It is seen from Eq.(9) that the necessary condition for formation 
of the fermion condensates is sufficiently strong four-fermion coupling g that gA/n 2 > 1, 
similar to the one in 4D NJL-model case [8]. In both cases, one can consider the four- 
fermion couplings as an effective interactions below some high momentum scale A. For 
the sake of satisfying the gap equation Eq. (8), the fine-tuning of the coupling constant 
g is also required. 

When T ^ 0, we must replace the vacuum expectation value SfcLi^^fc) by the 
thermal expectation value Y,k=i (^^a^t- m ^ ne real-time thermal field theory, this implies 
the substitution of the fermion propagator [9,10] 

-2n5(l 2 -m 2 ){l/+m) sin 2 9(1°, /1) (10) 



]/— m + ie 1/— m + ie 



with 



8111 U[L ' /ij = exp[/?(/° - /i)] + 1 + exp[/?(-Z° + W+ 1 ' ( j 

where m = m(T,fj) is the dynamical fermion mass at finite temperature T and finite 
chemical potential \x and the denotation (3 = 1/T. As a result, the gap equation at T 7^ 
becomes 

l = gl (12) 



dH 

(27T) 3 



I 2 — m 2 + ie 



-2n5(l 2 -m 2 )sm 2 6(l°,n) 



(13) 



A f, m A txT rT . T . . ,1 

= ^ I ^ arCtan m" 2A [ 2 ^'~ r)+ 2(l/ ' r)] J ( ^ 
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where we have used the denotations 



i 



o \/x 2 + y 2 exp^x 2 + =F r) + 1 



In 



1 + e 



with 



x — (5\ I |, y — f3m and r = 
Considering the gap equation (9) at T = 0, Eq.(12) may be changed into 



m(0) arctan — — = m arctan 1 — f(T) 

v ; m(0) m 2 JK J 



with 



f(T) =Tln{[l + e- (m ^ )/T 



1 + e 



-(m+n)/T 



]}■ 



(15) 
(16) 

(17) 
(18) 



It is easy to see that no matter m — fi > or m - /i < Owe always have f'(T) = 
df(T)/dT > 0, i.e. f(T) is a monotone increasing function of temperature T. To keep the 
left-handed side of Eq.(17) to be a constant, the first term in the right-handed side must 
decreases as T increases and finally has m — > in it at a critical temperature T c and a 
critical chemical potential fj, c . Hence taking m = in Eq.(17) we will obtain the equation 
to determine T c and ji c 



m(0) arctan 



A 



7T, 



m(0) 



ln(l + e^ c/Tc ) + ln(l + e '^ /Tc ) 



(19) 



Here the explicit momentum cut-off A is not important. In fact, as long as the ratio 
A/m(0) > 10 4 , then the value of arctan[A/m(0)] is almost equal to tt/2. Assume that is 
the case then the momentum cut-off A will disappear from the equation and we obtain 



m(0) = T c [ln(l + e^ c/Tc ) + ln(l + e '^ /Tc ) 



(20) 



We indicate that the /i — T criticality equation (20) comes from the two-point func- 
tion of frmions which corresponds to 0(1) order in 1/N expansion, hence no explicit 
Af-dependence is contained in it. It follows from Eq.(20) that 



T c /m(0) = 1/2 In 2, if /i c = 
/x c /m(0) = 1, if T c = 



(21) 



If scaled in m(0) by defining u c = /i c /m(0) and t c = T c /m(0) then the equation (20) of 
criticality curve of chemical potential and temperature can be rewritten as 



u c = 1 + 2t c ln[(l + VI - 4e-V*c)/2]. 



(22) 



These results obtained under the assumption A/m(0) ^> 1 are identical to the ones given 
by effective scalar field approach in the leading order of 1/N expansion [5] and also 
cosistent with the calculations of Ref. [4]. 

By means of the dynamical fermion mass, we may have more clear understanding of 
the feature of symmetry restoration phase transition. Assuming A/m(0) ^> 1 and in view 
of m < m(0), we will have arctan[A/m(0)] ~ arctan[A/m] ~ tt/2. Thus Eq.(17) will 
become A-indepedent, i.e. 



m(0) = m + T jln [l + e 



(m-n)/T 



+ In [l + e -( m+ ^ T ] } . 



(23) 
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When T ~ T c , we have m/T <C 1 thus Eq.(23) can be transformed into 



m(0) = T [in (l + e^ /T ) + In (l + e~^ T ) 



+ Tln 



1 + ( e m/r - l) + - 



In (l + e 



m/t^) + i n (i + e -M/T) 



+ 



(1 + e»l T ) (1 + e-^ T ) 

3\ 1 



2T 2 [1 + cosh(/x/T)] 



+ 



m 

J^3 



where we have made the power series expansion of [exp(m/T) — 1] in m/T and kept 
only the terms to the order 0(m 2 /T 2 ). Combining Eq.(24) with Eq.(20) (in which /z c is 
replaced by ji) we obtain 



m 



2r(l + cosh|) |T c ln 



2 1 + cosh 



Tin 



2 1 + cosh 



(25) 



m 



(24) 



By means of the approximate expression at T ~ T c 



In 



1 + cosh ^~ 
T 



In 
In 



l + cosh{ii + ^^ 

\t c t c t 



2 

2(1 + cosh -- 



// 



T, 



+ ^: T ^tanh " 



m (26) 



Eq.(25) may be changed into 



m 



2T (l + cosh 



In 



2 1 + cosh 



-^tanh^}(T c -T) when T Z T c 



(27) 

We see that the dynamical fermion mass m, as the only order parameter of the discrete 

symmetry breaking, will have (T c — T) 1 / 2 behavior at T ~ T c . This situation is similar 
to the one in 4 dimension NJL-model [10]. Hence, based on the same arguments given 
in Ref. [10], we can conclude that in the 3D Gross- Neveu model the discrete symmetry 
restoration at T > T c will also be a secon-order phase transition. 

In this model, no continuous symmetry breaking means no Goldstone boson appearing. 
An interesting question is that whether any (massive or massless) scalar bound state could 
exist. Such scalar bound state can be thought as the equivalent of a "Higgs" particle 
appearing in spontaneous breaking of a continuous symmetry. In present case, it could 
be the configuration SfcLi('^ fc V'fc)- However, we will prove that the answer to the above 
question is negative in 3D Gross- Neveu model. Based on the method taken in Ref. [11] in 
the real-time formalism of thermal field theory, we can obtain from the scalar four-point 
function of fermions the formal propagator for the "scalar bound state" expressed by 



r s (p) = -i/N(p 2 - Am 2 + is) 



K{p) + H(p) - iS(p) - 



R 2 (p) 



K{p) + H(p) + iS(p) 



(28) 



where the functions 



K(p) 



dx 



4tt 2 Jo [m 2 — p 2 x(l — x)] 1 / 2 
Am 2 + 4m y/jp + p 1 



arctan 



A 



[m 2 — p 2 x(l 



Zn^p 2 



In 



4m 2 — p 2 



, if arctan 



- x)] 1 / 2 
A 



7Y 



[m 2 — p 2 x(l — x)] 1 / 2 



- and < p 2 < Am 2 , 



5 



(2tt) 3 [[(/+p) 2 -m 2 ] 2 + e 
rf 3 l 

(27T) 



5(p) = 2tt 2 / -^-^5(Z 2 - m 2 )5[(Z + p) 2 - m 2 ] [sin 2 9(1°, fi) cos 2 0(Z° + p°, //) 



+ cos 2 #(Z°,/i) sin 2 ^(/°+p°,/i) 

J R(p)= 7r 2 / -^ ( 5(/ 2 -m 2 ) ( 5[(/+p) 2 -m 2 ]sin2^(/ , At )sin2^(/ +p ,/i). (29) 
7 (27T)" 1 

It is indicated that when deriving Eq. (28) we have used Eq. (12). As a result, the A- 
dependent divergent sectors in T s (p) have cancelled each other and this is just a peculiarity 
in the leading order of 1/N expansion. Formally p 2 = 4m 2 seems to be a mere simple 
pole of T s (p), but this is not true. First of all we notice that when p 2 — > 4m 2 , 

lim K(p) ^oo (30) 
and K{p) has a logarithmic singularity. Next, we can write 

„, x 1 f d 2 l p 2 + 2u l p°-2\ T || P |cosy? p 2 -2c<V ) + 2| / || P |cosy? N 

= 3-2 / — U— — :+— — -) 



871 J [p2 + 2uj lp o _ 2| 1 || p \cos(f} 2 + e 2 [p 2 - 2^ + 2| I \\ P | cosy?] 2 + e 2 ' 

1 ^2 1/2 

- + (p°^-p ,/^ -//)}, "i = (Z +m 2 ) . (31) 



exp[/3(cj/ -/x)] + 

The angular integrations in Eq. (31) are of the following general forms and results: 

/ d ^TZ 4-^ = °' ( 32 ) 

Jo o + ccos(p±e z 

where b and c are constants. Hence the function H(p) = 0. As for S(p) and R(p), when 
p 2 = Am 2 each of them may be generally expressed by 

A(p)\ p2=4m2 = Jd 3 l5(l 2 -m 2 )5[(l+p) 2 -m 2 ]f(l ,p ^)\ p2=4m2 

r°° d\ 1 || 

'0 Aooi 



r dl * 11 1 1 F dtp 5(2m 2 - cu lP ° - I 7 || P I cos^)/(-^,p ,/i), 

JO 4W; JO 



(33) 

where f(l°,p°,fi) is a finite function. In Eq.(33), the condition |cos</?| < 1 limits | / 

to the only possible value | I \ — \ P |/2. Let tan(</?/2) = t, then the integration of the 
angle ip in Eq. (33) can be written by 

J + J J dip 5(2m 2 - uip° - I I || V I cos <p) 
= 2^°°rft |5[(2m 2 - uj iP ° + I T || P \)t 2 + (2m 2 - uj iP ° - | 7 || P |)] 

+5 [(2m 2 - - I 7 || P |)t 2 + (2m 2 - uj iP ° + | 7 || P |)]} . (34) 



6 



Since the (^-functions have non-zero values only if | I \ — \ P |/2, we can set 

1 



and obtain 



which further lead to 



I | = (-±e)|p|, e = + 



2m 2 - uip° + | / || P 



e 2 P 



— -±2 



2m 2 - ujip - | / || P 



p +4m 2 

^2 

-P 



(35) 



(36) 



B = 2]im£- - { I dt8 

s^O I I |,|P|/2 I 7 



g 1 - P 



P +4m 2 



+ / dt 5 
Jo 



= 2 (^)^ im W =2 (^ )1/2 lim 



^2 



£^0 



2^ 4 

e 2 P 



P +4m 2 

« l,|P|/2 



e P 



P 



UHIpI/2 



/ I - I P 1/2 



(P +4m 2 ) 1 / 2 - 
= } 5{\ I | - | P |/2). 

P 

Substituting Eq. (37) into Eq.(33) we will have 



(37) 



Mp)\ p 



,2_4 m 2 



21 P 



! /(-Vi 2 + m 2 ,p ,/i) 



(38) 



|i|=|P|/2 



This means that both S{p)\ p 2 =Am 2 and R(p)\ p 2=4 m 2 with the form of Eq. (38) are finite. 
To summarize the above results, we can obtain 



lim (p 2 — Am 2 )Ts(p) = 0, 

p 2 — +4m 2 



(39) 



hence p 2 = 4m 2 is not a simple pole of Ts{p) and Ts{p) does not represent the propagator 
for any scalar bound state. This proves that no scalar bound state could exist in this 
model. This result, seen from Eqs. (28) and (30), is true at T = 0, and also maintained 
at a finite temperature. We point out that the above result at T = is different from the 
one obtained in Reference [5]. 

In conclusion, in 3 dimension Gross-Neveu model, the fermion codensates induced by 
the scalar four- fermion interactions will spontaneously break the discrete Z 2 , the spe- 
cial parity V\, V2 and the time-revesal symmetries and lead to generation of dynamical 
fermion mass. By means of the real-time thermal field theory in the fermionic bubble 
diagram approximation, it is clearly shown that the dynamical fermion mass m, as the 
order parameter of symmetry breakings, will decrease as temperature increases and finally 
becomes zero at some critical temperature T c . This implies restorations of these discrete 
symmetries at T > T c . The equation of the criticality curve of chemical potential and 
temperature will be indepedent of the momemtum cut-off and identical to one obtained 
by the auxialiary scalar field approach in the 1/N expansion if the ratio between the 
momentum cut-off and the dynamical fermion mass at T = is sufficiently large. Similar 
to the case in 4D NJL-model, the dynamical mass m has the (T c — T) 1 / 2 behavior at 

T ~ T c and this fact indicates the second-order phase transition feature of the symmetry 
restorations. We have also proven by explicit analysis of the fermionic scalar four-point 
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functions that there is no scalar bound state in this model. 
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